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We study a limit of the nearly-Peccei-Quinn-symmetric Next-to-Minimal Supersymmetric Stan- 
dard Model possessing novel Higgs and dark matter (DM) properties. In this scenario, there nat- 
urally co-exist three light singlet-like particles: a scalar, a pseudoscalar, and a singlino-like DM 
candidate, all with masses of order 0.1-10 GeV. The decay of a Standard Model- like Higgs boson to 
pairs of the light scalars or pseudoscalars is generically suppressed, avoiding constraints from collider 
searches for these channels. For a certain parameter window annihilation into the light pseudoscalar 
and exchange of the light scalar with nucleons allow the singlino to achieve the correct relic density 
and a large direct detection cross section consistent with the CoGeNT and DAMA/LIBRA preferred 
region simultaneously. This parameter space is consistent with experimental constraints from LEP, 
the Tevatron, T- and flavor physics. 



The Next-to-Minimal Supersymmetric Standard 
Model (NMSSM) is a well-motivated extension of the 
Minimal Supersymmetric Standard Model (MSSM) by 
a gauge-singlet chiral superfield N, designed to solve 
the /i-problem of the MSSM. Its superpotential and soft 
supersymmetry-breaking terms in the Higgs sector are 

W = ANH u H d + ^N 3 , 



Vsoft = m 2 Hd \H d \ 



m 2 Hu \H u \ 2 



■m 2 N \Nf 



- (\A\H u H d N + h.c.) + {^A K N S + h.c.) . (1) 

Here H dl H u and N denote the neutral Higgs bosons 
corresponding to Hd, H u and N, respectively. 

In this work, we examine an NMSSM limit given by 
two conditions. The first one is n <C A which is protected 
by an approximate Peccei-Quinn (PQ) symmetry. It is 
well-known that a light pseudoscalar a\ will be generated 
by the spontaneous breaking of such a U(l) symmetry 
(the phenomenology on a light a\ has been thoroughly 
studied in the R-symmetry limit [U[2]). As noted in [3], 
at tree level the PQ limit implies an upper bound on the 
lightest scalar mass approximately proportional to 
A 2 . Here we address the further limit of A < 0.1, lead- 
ing to the simultaneous emergence of a light singlet-like 
scalar hi and a light singlino-like lightest superpartner 
Xi- For mildly small values of A (A > 0.05) studied in 
this letter, typically A(Aqut) ~ 0(0.1), a natural order 
for a perturbative parameter. We stress that this sce- 
nario differs from the light a\ case of [TJ [2], in that hi, 
ai, and Xi are a ll °f order 0.1 — 10 GeV. It also differs 
in that decays of the Standard Model (SM)-like Higgs 
boson to hihi and aiai pairs are generically suppressed. 
Thus hi and ai are hidden from four-fermion searches at 
LEP [4 and the Tevatron [5] designed to test a light a\ 
scenario. Meanwhile, due to annihilation into ai and ex- 
change of hi , for a certain window of the parameters, the 



correct relic density and a large spin-independent (SI) di- 
rect detection cross section consistent with the CoGeNT 
and DAM A/LIBRA preferred region can be achieved for 
the DM candidate xi- Therefore, we refer to this limit 
as the "Dark Light Higgs" (DLH) scenario. 

We begin with an analysis of the light spectrum in the 
DLH scenario. For convenience we define two parameters 



A/i , , A x 

mz H tan p 



(2) 



with ji = A (TV), e has an impact on Higgs physics, as 
exhibited in FIG. [I] In the first column of FIG. [I] we plot 
m h 1 ,a 1 ai a g ams t £ for a random scan as defined in the 
figure caption. NMSSMTools 2.3.1 and MicrOMEGAS 
2.4.Q [6j[7] ar e our analysis tools used in this letter. 

The scan results in FIG. [I] can be understood analyti- 
cally as follows. Because of the spontaneous breaking of 
the approximate PQ symmetry, a\ is a pseudo-Goldstone 
boson and its small mass m 2 ai « — 3kA k /jl/X is protected. 
For Xi, k <C A <C 1 implies that it is dominantly singlino 
and its mass is m Xl « v 2 X 2 sin 2/3/ ja + 2/^/i/A, where 
v = 174 GeV and tan/3 = (H u )/(H d ). For A < 0.1, 
fx of order a few hundred GeV, and «/A on the order of 
a few percent, m Xl drops below 10 GeV. 

More interesting is the CP-even spectrum. For analytic 
convenience we consider moderate tan /?, although the 
qualitative properties of the figures are also present for 
lower tan (3. In the small A + PQ limit hi has a mass 
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at tree level. The heaviest state is strongly down-type, 



with a mass mi 

fl3 



m^ Hd ~ A\ (where the minimization 



condition for v d is used), and the middle state is SM-like. 

The hi mass is also lifted by quantum corrections, and 
the strong singlet-like nature of hi suppresses contribu- 
tions from all particles running in the loop except Higgs 
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FIG. 1: Masses of hi (top-left), a\ (middle-left), and %i 
(bottom-left); branching ratios of h<2 into ftifti (top-right) 
and a±ai (middle-right), and correlation between e and e' 
(bottom-right). Points are taken randomly from the ranges 
5 < tan/3 < 50, 0.05 < A < 0.5, 0.0005 < k < 0.05, -0.8 < 
e' < 0.8, -40 GeV < A K < 0, and 0.1 TeV < fi < 1 TeV. 
(As an illustration, we assume soft squark masses of 1 TeV, 
slepton masses of 200 GeV, A u ^, e parameters of 750 GeV, 
and bino, wino and gluino masses of 100, 200 and 660 GeV, 
respectively, for all numerical analyses in this letter.) Green 
points cover the whole scan range, red points correspond to 
A < 0.30, k/X < 0.05 and /x < 400 GeV, and blue points 
correspond to A < 0.15, k/X < 0.03 and ji < 250 GeV. 

bosons and Higgsinos. Setting e — » for these loop dia- 
grams, we find an uplifted singlet mass in the MS scheme 



2 AV /i 2 tan/3 3 



(4) 



Fixing all other parameters, the upper bound on m\ is 
achieved for e —> and is lowered to about or below 10 
GeV in the small A + PQ limit. 

On the other hand, increasing e rapidly decreases • 
Vacuum stability (i.e., (m| ) tre e + ^ m h 1 > 0) indicates 
an upper bound on e 2 
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In the small A + PQ limit and for natural values of /i, 
kmo^l is small. This fact will be relevant for collider 
constraints discussed below. The right-bottom panel of 
FIG. [I] also shows that A\ is usually close to /itan/3 for 
blue points, so we will take a smaller range of e' in our 
DM analysis. 

The tree-level mixing parameters of the light scalar are 

Su«-^U + ^), S luf »^, (6) 
/i tan p \ rriz ) thz 



indicating a mostly singlet /down admixture in the limit 
e —> and an approximately pure singlet (i.e., Si s — » 1) 
in the further limit of small A or large tan/3. 

There are three main processes by which present ex- 
periments potentially constrain this scenario: (1) decays 
of the SM-like Higgs h 2 through ftifti and a\a\, (2) T(ns) 
decays to jhi or jai , and (3) flavor physics. 
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FIG. 2: Constraints from the decays ft 2 —> ftifti — >• 4/ (top) 
and from the decays T — >> 7fti(fti -4- /jl/jLjIttt, KK) (bottom). 

= <7h, 2 Br(ft2 — >• ftifti — >• 4/i). To show the constraint 
from the 2/i2r channel on the same plot we convert it into 
an effective constraint on 4/x by rescaling it with Br ^ 1 



Br(/ii — >tt) 

(a model-independent quantity). A^ is a tree-level coupling 

of the down- type interaction *y^ d hifdfd- Gray and blue 

points correspond to the gray and blue points in FIG.[3|Pur- 
ple bands correspond to the points in the scan of FIG~]4] 

Similarly to the light a\ scenario of [2 , relevant con- 
straints may come from the searches for [4j |5] 

ft 2 -> ftifti, aiai -> 46, 4r, 2b2r (LEP), 
ft 2 — > ftifti, aiai -4- 4/i, 2/i2r (Tevatron). 

However, in our case the tree-level couplings of ft-2 to ftifti 
and aiai are suppressed. This can be seen as follows. 
Since fti is strongly singlet-like and is up-type, the 
coupling Uh 2 h 1 h 1 is (f° r a complete formula, see [TO] ) 
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Here we use the mixing parameters at lowest order in e 
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for moderate tan/3. Similarly, one can find yn 2 a 1 a 1 = 
Vh 2 h 1 h 1 at this order. Both Br(ft2 —> ftifti) and Br(ft,2 — > 
a\ai) are thus suppressed by Xe <C 1, as is shown in the 
right column of Fig. [I] (Instead, h 2 can dominantly de- 
cay into xi an d X2, while %2 dominantly decays into light 
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Higgs bosons and Xi- These facts imply rich Higgs phe- 
nomenology in the DLH scenario and can dramatically 
change the strategies of searching for the SM-like and 
light Higgs bosons at colliders [11.) The asymmetry in 
Br(/i2 — >• h\h\) w.r.t. e is caused by an 0{e 2 ) correction 
with the opposite sign of the term in Eq. ([7]). 

The Tevatron constraints from the search for hi — >• 
h\h\ — >• 4/ are illustrated in the upper panel of FIG. |2j 
Almost all points survive. Similar limits from LEP are 
avoided easily for the present parameter values, because 
777^2 is above the kinematic threshold 1 . 

T physics constrains models with light states through 
T — >• 7(/ii,ai) —> 7(MM, tttt, KK). Fig. [2] shows the con- 
straints from searches for these decays on the effective 
coupling Ad of the light state to down- type fermions [8j|9]. 
At tree level, « ^(A + f^), and the scan points typ- 
ically approach the constrained region only for A > 0.15. 

5-physics may also add non-trivial constraints with a 
light a\ (e.g., see [10]) or h±, because flavor- violating ver- 
tices 6(d, s)(ai, hi) can be generated at loop level. These 
vertices, however, depend strongly on the structure of 
soft breaking parameters (e.g., see [H]). For the input 
parameters to NMSSMTools used in the scan, the points 
in the figures are consistent with all 5-physics constraints 
including B s — » /i/i, B& — » X s \i\i^ b — >• 57, etc. In addi- 
tion, though not included in NMSSMTools, we also check 
the constraints from D meson decays (e.g., D — »• / + /~). 
Because of the singlet-like nature of hi and ai, D-physics 
constraints are very weak and can be satisfied easily. 

To study the DM physics in the DLH scenario, we 
perform a second random scan over its parameter re- 
gion (a narrower region than the one in the first scan). 
FIG. [3] shows that the xi DM candidate is character- 
ized by a larger spin-independent direct-detection cross 
section asi, compared with typical supersymmetric sce- 
narios. For certain parameter window, the correct relic 
density and a large <jsi consistent with the CoGeNT and 
DAM A/LIBRA preferred region [14] can be simultane- 
ously achieved, and the scenario remains consistent with 
current experimental bounds (particularly from flavor 
physics and Higgs searches). This has been considered 
difficult or impossible in supersymmetric models [T5HT7] . 

The large <jsi is mainly due to the h\— mediated 
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FIG. 3: Cross section of SI direct detection for xi- The 
scan is over all parameters, in the ranges 0.05 < A < 0.15, 
0.001 < k < 0.005, \e'\ < 0.25, -40 < A K < GeV, 
5 < tan/3 < 50 and 100 < fi < 250 GeV. The dark blue 
(dark) points have a relic density 0.09 < Qh 2 < 0.13. The 
red contour is the CoGeNT favored region presented in |13| 
and the two blue circles are the most recent interpretations 
of fitting CoGeNT + DAM A/LIBRA [14]. All contours as- 
sume a local density which may be sensitive to the relic den- 
sity. The purple, brown, and black lines are the limits from 
CDMS [18], CoGeNT [13], and XENON100 [19], respectively. 
Most CoGeNT favored regions have a tension with the CDMS 
constraints. Consistency between the CoGeNT preferred re- 
gions and the XENON100 constraints can be achieved within 
the scintillation-efficiency uncertainties of liquid xenon [14] . 



t— channel scattering xiQ ~^ X1Q1 an d 0"si 
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The hiXiXi coupling is reduced to y^xixi ~ — for 
a singlino-like xi an d singlet-like hi. The dependence of 
a si on m^ 4 is illustrated in the left panels of FIG. [4] For 
the parameter values given in the caption, the LEP search 
for hi — >• bb sets the lower boundary of the contoured 
region, flavor constraints control the upper-right, vacuum 
stability sets the upper- left limit, and the upper bound 
on the relic density controls the left and right limits. The 
sensitivity to tan f3 enters mainly via • 

The xi renc density is largely controlled by the ai- 
mediated annihilation X1X1 ~~ ^ //? with cross section 



1 The LEP and Tevatron constraints from the channel /12 — > a\a\ 
are included in NMSSMTools and in our code, respectively. 
Points are omitted if the limit is violated. Similarly, the con- 
straint from T — y 7<ji is checked by NMSSMTools, so we present 
only the limit from T — > 7/11 in FIG. [2] For the numerical 
results presented in this letter we incorporate all built-in checks 
in NMSSMTools 2.3.1 (including those from LEP Higgs searches, 
superpartner searches, — flavor physics, Z-decay, 775 physics, 
etc.), exc ept t he DM relic density. The difference between FIG.[l] 
and FIG. |3|4|is that in the latter, Qh 2 < 0.13 is also required. 



3| y aiXlXl y aif f\ 2 (l - m 2 f /m 2 Xl )^ 2 
°ff v xi ~ 7 9 a , (10) 
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v Xl denoting the relative velocity of the two Xi s 0- $v xi -^0 
reflects the deviation of 2m Xl from the ai resonance. In 
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FIG. 4: Contours of asi (top-left), Qh (top-right), rrih 1 (left- 
bottom) and 5 Vxi ^o (right-bottom) on the \i — tan/3 plane, 
with A = 0.12, k = 2.7 x 10" 3 , e' = 0.15 and A K = -24 GeV. 

the typical case m ai > 2m Xl > Im^, the relic density is 



Qh 2 



n -| / m ai \ ( r ai \ / 0.003 A 2 /ai m n 

u,J -Vl5GeV/ ^lO-SGeVy V^ixixi/ V A v y 

(Sfv / ^^^)/erfc(2.2) 
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where Xf = m Xl /Tf is the freeze-out point. As a measure 
of thermal suppression, 5 Vxi ->o enters the complementary 
error function obtained from the integral over the Boltz- 
mann distribution. The inverse dependence of Qh 2 on 



is shown in the right panels of FIG 



Its sensi- 
tivity to fi is mainly through 5 Vxi ^o, as m x Jm ai ex ^JJl 
for tan j3 > 5. To achieve the correct relic density requires 
5 Vxi ->o « 0.30 — 0.35, which implies A K & — 3.5m Xl , with 
a tuning range about ±0.1m Xl . We emphasize that this 
process does not generate an antiproton or 7-ray flux 
in tension with existing cosmic-ray data because of the 
Breit-Wigner suppression effect today [20] . 

Finally, a benchmark point corresponding to the stars 
in FIG. and FIG. [3] is given in Table [TJ We would like to 
point out that the chosen set of parameters in the squark, 
slepton and gaugino sectors in this letter provide a real- 
ization of the DLH scenario. Changing them can change 
the details of the phenomenology, but the basic features 
will remain intact. We reserve an extended phenomeno- 
logical analysis of this scenario for future work. 
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TABLE I: Benchmark point. We use the units cm 2 for <tsi and 
GeV for dimensionful input parameters, and denote Br(/i 2 — »• 
hihi) as Brhh and Br(/i2 —> a±ai) as Braa. Soft sfermion and 
gaugino parameters are as given in the caption of FIG.[l] 
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